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Abstract 

The concept of weak ergodicity breaking is defined and studied in the context of deterministic 
dynamics. We show that weak ergodicity breaking describes a weakly chaotic dynamical system: a 
nonlinear map which generates subdiffusion deterministically. In the non-ergodic phase non-trivial 
distribution of the fraction of occupation times is obtained. The visitation fraction remains uniform 
even in the non-ergodic phase. In this sense the non-ergodicity is quantified, leading to a statistical 
mechanical description of the system even though it is not ergodic. 
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INTRODUCTION 



The relation between deterministic non-linear dynamics and statistical mechanics is a 
topic of wide interest l|, l2|. In particular, the domain of validity of the ergodic hypothesis 
for low dimensional dynamical systems is fundamental in our understanding of the foun- 
dations of statistical mechanics. In usual setting of statistical mechanics the phase space 
of a closed system is divided into equally sized cells, then if the system is ergodic it will 
spend equal amount of time in each cell in the limit of long measurement time. Ofcourse 
many dynamical systems are non-ergodic, for example the phase space may decompose into 
regions, where the system cannot access one part of the phase space if it started in a different 
part. Such non-ergodicity is called here strong non-ergodicity. 

A second type of non-ergodicity, called weak non-ergodicity is investigated in this 
manuscript. As conjectured by Bouchaud jsl in the context of glassy dynamics, in weak 
ergodicity breaking the phase space is a priori not broken into mutually inaccessible re- 
gions. However the system does not spend equal amount of time in equally sized regions 
of phase space (see detailed definition below). While occupation times of cells in phase 
space remain random, we showed in the context of stochastic dynamics that distribution of 
occupation times is universal {4] in a way that will be made clearer later. The goal of this 
paper is to investigate weak non-ergodicity for a deterministic process. 

In terms of stochastic theories like the trap model, the random energy model Sj and the 
continuous time random walk [4] it is known that when the dynamics exhibits aging and 
anomalous diffusion it may exhibit weak ergodicity breaking. It is also well known that a 
wide variety of deterministic dynamical systems exhibit anomalous diffusion 

mm . Thus 

the question arises: Can deterministic dynamical systems exhibit weak ergodicity breaking? 
In this manuscript a detailed definition of weak ergodicity breaking (WEB) in the context 
of deterministic dynamics is introduced. It is shown that that a low dimensional model, a 
deterministic non linear map which has no disorder built into it, yields WEB. Non-trivial 
(in a way defined later) probability density function (PDF) of fraction of occupation time 
is shown to describe equilibrium properties of the dynamical system. In this sense the 
non-ergodicity is quantified, enabhng us to perform statistical mechanical description of the 
system although it is not ergodic. 
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GEISEL'S MAPS 



Probably the simplest deterministic maps which lead to normal and anomalous diffusion 
are one dimensional maps of the form 

xt+i = xt + F (xt) , (1) 

with the following symmetry properties of F (x) : (i) F (x) is periodic with a period interval 
set to 1, i.e., F (x) = F {x + N) , where N is an integer, (ii) F (x) is anti-symmetric; namely, 
F (x) = —F {—x) , note that t in Eq. ([T]) is a discrete time. Geisel and Thomae [9| considered 
a wide family of such maps which behave like 

F (x) = ax' for x 0+ (2) 

where z > 1. Eq. ([2]) defines the property of the map close to its unstable fixed point. In 
numerical experiments introduced later we will use the following map 

F(x) = (2x)^ 0<x<l/2 (3) 

which together with the symmetry properties of the map define the mapping for all x. The 
functional form of the map on a unit interval was introduced by Pomeau and Manneville 

1^. Variations of these maps have been investigated by taking 



to describe intermittency 



into account time dependent noise jn], quenched disorder [l^, and additional uniform bias 



which breaks the symmetry of the map [13|- We consider the map with L unit cells and 
periodic boundary conditions. In Fig. [1] we show the map, and one realization of a path 
along three unit cells. 

The dynamics of the map in the vicinity of the unstable fixed points, is governed by a 
q exponential sensitivity of the trajectories to initial conditions, i.e., weak chaos. The role 



of such q exponential in mode 
a subject of ongoing research 



s re. 



ated to ours (the logistic map) and Tsallis statistics is 
14, Q, Q. ideas related to 

Tsallis statistics are wrong due to misconceptions on basic notions of ergodic theory, here we 
show that for the weakly chaotic map under consideration, the correct statistical framework 
is given by the concept of WEB. 

In an ongoing process a walker following the iteration rules may get stuck close to the 
vicinity of one of the unstable fixed points of the map. It has been shown both analytically 



3 




FIG. 1: The mapping Eqs. (jll3p with z=3 and L=3, the dashed hne shows the curve xt+i = xt- 
Fixed points are at xt = 0, 1, 2. A specific reahzation of the dynamics is also drawn (thin line). 

and numerically, that the PDF of escape times from the vicinity of the fixed points decays as 
a power law [9]. To see this, one considers the dynamics in a half unit cell, say < x < 1/2. 
Assume that at time t = the particle is on x* residing in the vicinity of the fixed point 
X = 0. Close to the fixed point we may approximate the map Eq. ([1]) with the differential 
equation dx/dt = F (x) = ax^ . This equation is reminiscent of the equation defining the q 
generalized Lyaponov exponent [l6]. The solution is written in terms of the q exponential 
function, exp^ (y) = [1 + (1 — g) y]^^'-^"''^ where q = z and 

X (t) = X* exp^ \ax*^^-^h] . (4) 



We invert Eq. (j4j) to obtain the escape time from x* to a boundary on b{x* < b < 1/2) which 
b ^ 

is t ~ / [F{x)]~^dx, using Eq. we obtain t ~ (1/a) [{x*^-^)/{z - 1) - {b'^-')/{z- 1)] a 

X* 

luq ( i.e., the inverse function of exp^ ) behaviour. The PDF of escape times ip (t) is related to 
the unknown PDF of injection points rj (x*) , through the chain rule (t) = rj (x*) \dx*/dt\ . 
Expanding r] (x*) around the unstable fixed point x* = one finds that for large escape 
times 

V' w ~ pfTT^r a = iz-ir' (5) 

|r(-a)| 

where A depends on the PDF of injection points, namely on how trajectories are injected 
from one cell to the other. The parameter A is sensitive to the detailed shape of the map 
while the parameter z depends only on the behaviour of the map close to the unstable fixed 
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joints. When z > 2 the average escape time diverges [5|, l9|, ll8| , and in turn yields aging 



191], non-stationarity 



20|, anomalous diffusion 



21 



2211 and as we show here WEB. 



WEAK ERGODICITY BREAKING 

The map is said to be ergodic if for almost any initial condition (excluding paths starting 
at fixed points) the path spends equal amounts of time in each cell in the limit of long 
measurement time. Namely jt = 1/L, where is the time spent in one of the lattice 
cells and t is the total measurement time. Another method to quantify the dynamics is 
to consider the fraction of number of visits per cell n^/n where is the number of visits 
in one of the lattice cells and n is the total number of visits (or total number of jumps 

oo 

between cells). For maps with finite average escape time (r) = J rip (r) dr, i.e., z > 2, the 



total measurement time t ^ n (t) and the occupation time of one cell ~ (r) hence 
n^/n = 1/L, the equalities above should be understood in statistical sense. We call n^/n 
the visitation fraction, which plays an important role in WEB. 

More generally the fraction of occupation times (visitation fraction) in m < L cells is 
t"^/t {n"^/n), where t"^ (ra'") is the occupation time (number of visits) in m cells respectively. 
For ergodic systems the PDF of the fraction of occupation time is 

fC-).,C--'-), (6) 

and the visitation fraction is 



t J \ t L 



in statistical sense, both equalities are valid only in the limit of long measurement time 
(t — > oo). Note that the specific choice of the m cells is not important. Maps exhibiting 
strong non-ergodicity do not obey Eqs. (|6|7|) . Dynamics obeying Eq. ([7]) but not Eq. ([6]) 
are said to be weakly non-ergodic. As we will show when z > 2 WEB is found. 

We start our analysis by discussing numerical simulations of the map with system of size 
L = 9. In Fig. [2] (a) we show two paths generated according to the mapping rule Eqs. 
with z=3. As one can see each path gets stuck at one of the cells for a time which is of the 
order of the measurement time. During the measurement the path visits all the L cells many 
times. In Fig. [2] (b) we consider the case z = 1.5, we show that the fraction of occupation 
time obeys Eq. ([6]) with m = 1 namely t^/t = 1/L, with L = 9. On the other hand in Fig. [2] 
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FIG. 2: (a) Two paths (dashed and sohd curves) generated according to the mapping rule with 
z=3, for each path the particle gets stuck in one of the cells for a time which is of the order of the 
measurement time, (b) Histogram of the fraction of occupation for a path with z=1.5. For this 
ergodic case the fraction of occupation times in the cells are equal. (c)+(d) The same as in (b), 
but this time for the non-ergodic phase with z=3. 

(c,d) we show the fraction of occupation time histogram for two initial conditions obtained 
using z = 3. Since each path is localized in one of the cells for a time which is of the order 
of the total measurement time, the histogram for this case looks very different than that of 
the ergodic case and clearly Eq. ([6]) is not valid. Histograms of the visitation fraction in 
each cell are presented in Fig. [3l for both ergodic {z = 1.2) (a) and non-ergodic {z = 3) (b) 
maps. The Fig. demonstrates that for both cases Eq. ([7]) holds, i.e., the visitation fraction 
in each cell is given by 1/L. From Fig. [2] we see that the fraction of occupation time is 
random, hence we investigate its distribution. We will find universal PDFs of occupation 
time, which generalize the ergodic rule Eq. 

The PDF for t™( the occupation time of m cells) in the case of diverging mean escape 
time is now derived. In particular we show that the assumption Eq. ([7]) with power law 
waiting times, leads to nontrivial PDFs of fraction of occupation times of the map. We 
denote by /°"* (t^) the PDF of t"*, in the case that the path isn't within the m cells at time 
t and that during the measurement time t there were n jumps. Then 

fZt (n = {t"' - E ^r) / (T„ < t < T„+o) , (8) 

where the summation in the delta function is over all sojourn times within the m cells, 
namely is the k'th sojourn time within the m cells, T„ is the time after n jumps and / {x) 
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FIG. 3: Histogram of the visitation fraction in each of the 9 lattice cells, (a) z=1.2 dashed and 
solid lines correspond to two initial conditions, (b) The same as (a) but now z=3. In both cases 
the visitation fractions in each cell are equal, with small fluctuations. 

is equal 1 if x is true and otherwise, the brackets denote average over all sojourn times r. 
Double Laplace transform (LT) of Eq. ([8]) yields /°"* (u) = 

e-^'^^e-^'S (u - i: r™) J (T„ < t < T„+i) dUdt) 

^0 \ k=l J I 

where ip (s) is the LT of {t). We assumed that waiting times in cells are not correlated. 
This and other assumptions are checked later using numerical simulations. In a way similar 
to Eq. dHD we denote by (t™) the PDF of t"^ in the case that the particle is within the 
m cells at the end of the measurement, then in double Laplace space we find 

fZ (u) = {u + s) r-- is) ^-r^^'Ju + s) ^ 

u + s 

Considering an ensemble of initial conditions uniformly distributed in a cell, the probability 
for the path to be within the m cells at the end of the measurement is m/L, thus the double 
Laplace transformed PDF of the occupation time of the cell, given that there were n visits 
during the measurement is 

fn,s (U) = jfZ (U) fZ (U) . (11) 
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Using Eqs. (pfTOl) and the assumption /n = m/L Eq. ©we rewrite fn,s (u) as 

fn,s (u) 



ml — ip {u + s) ^ L — ml — ip {s) 



L s + u 

+ (12) 

Let ft (t™) be the PDF of occupation time of the m cells, its double LT is given by 

oo 

fs (u) = ^ fn,s 



n=0 



ml — ip {u + s) ^ L 



ml — ip (s) 



s + u 



" L — 771 , > 



1 _ ^mjL _^ _g) ^ 

Taking the limit m, s ^ (using the asymptotic behaviour of the escape time PDF ^ip (s 
1 — As'^^\ and inverting the double LT one obtains Lamperti PDF 
occupation time 



(13) 



23( 1 of fraction of 



sm 



(14) 



where R 



(1 - ) + (*^) + 2i? (i - ^) (^) cos 

which is valid only when t is large. Eq. (|T4l) for m = 1 is in agreement 
with the formalism developed in ^. Here we showed that the WEB condition n"'' /n = m/L 
(i.e., the visitation fraction is equal to the probability of occupying the m cells in ensemble 
sense) and the power law behaviour of the sojourn time PDF lead to the statistical law Eq. 

The assumptions made in the derivation of Eq. (IT^ are tested using numerical simulations 
of the map. Plots of the PDF of fraction of occupation time for three different values of 
z are shown in Fig. HI We used 10^ trajectories and measured the occupation time of 
m = 3 cells in a system of size L = 9. In the ergodic case z < 2 the PDF is just a delta 
function centered at m/L (see Fig. H] (a)). In the non-ergodic cases the maxima of / {t"^/t) 
are at and 1. These events t"^/t ^ l{t"^/t ^ 0) correspond to trajectories where the 
particle occupies (does not occupy) one of the observed cells for the whole duration of the 
measurement, respectively. In the non ergodic phase z > 2 two types of behaviors are found. 
For z = 2.25 the PDF of fraction of occupation time has a W shape, with a peak in the 
vicinity of m/L. While for 2; = 3, we find a U shape PDF and the probability of obtaining 
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FIG. 4: (a) The PDF of the fraction of occupation time of m=3 cells out of L=9 cells for z=1.5, 
In this ergodic case the PDF is very narrow and centered around m/L. (b) Same as (a), however 
here z=3. The PDF has a U shape, with maxima at and 1. (c) Same as (a) with z=2.25. In this 
case a W shaped PDF is obtained, with a peak at the vicinity of m/L, but yet the maxima are at 
and 1. The curve is Eq. (jl4p without fitting. Only the parameter z determines the shape of the 
PDFs, while detailed shape of the map is unimportant. 

fraction of occupation time equal to m/L is almost zero. Comparison of theory Eq. iHM 
with the numerical simulations of the map shows excellent agreement without any fitting, 
thus justifying our assumptions. Note that if we limit our model to integer z (i.e., assuming 
an analytical behavior) then we either observe WEB {z > 2) or the special border point 
between WEB and ergodicity z = 2. The latter point exhibits very slow convergence and 
logarithmic corrections which should be the subject of future work. 

CONCLUSION 

A detailed characterization of weak ergodicity breaking was established as follows: (i) 
The phase space is not divided into mutually inaccessible regions, thus the system visits 
all cells in its phase space for almost all (all but a set of measure zero) initial conditions, 
(ii) The visitation fraction fluctuates very slightly in the limit t ^ oo (note that n and ra™ 
exhibit strong fluctuations), (iii) The visitation fraction is determined by the probability of 
finding a member of an ensemble of particles within the given cell, (iv) The total time the 
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system spends in each cell fluctuates strongly, and is described by Eq. (fT^ . implying that 
statistics of occupation times depend only on the parameter z describing the non-linearity 
and not on the shape of the map far from the fixed points (e.g. a in Eq. ([2])). A relation 
between weak chaos, deterministic anomalous diffusion and weak ergodicity breaking was 
found. Since the former behaviors are wide spread in dynamical systems, it is possible that 
the non-ergodic dynamical scenario might emerge in other fractional dynamical systems, 
and then the present concepts naturally lead to a new type of weakly non-ergodic statistical 
mechanics. 
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